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Abstract 

A simple but empirically efficient heuristic algorithm for the edge-coloring of graphs is 
presented. Its basic idea is the displacement of "conflicts" (repeated colors in the edges 
incident to a vertex) along paths of adjacent vertices whose incident edges are recolored 
by swapping alternating colors (that is, doing a Kempe interchange). The results of 
performance tests on random cubic and A-regular graphs are presented, and a full 
implementation of the algorithm is given to facilitate its use and the reproducibility 
of results. 

1 Introduction 

Deciding if a graph can be A-edge-colored (that is, whether or not it is of class one) is 
known to be A^P-complete even for maximum degree A = 3 (Holyer [5J). Therefore, assu- 
ming P 7^ NP, the problem of 3-edge-coloring a cubic graph becomes intractable as the 
number of its vertices increases. Possible ways to "break" intractability are restricting the 
problem or admitting probabilistic solutions. Here we take the second approach, proposing 
a simple and fast heuristic algorithm to 3-edge-color a cubic graph G. Experimental data 
suggest that, when G is 3-edge-color able (an asymptotically almost sure property by the 
results of Robinson and Wormald [8]), our algorithm efficiently gives an explicit solution 
and, hence, allows us to break intractability. The algorithm is easily generalized to graphs 
with maximum degree A. All graphs are assumed to be connected and with no loops or 
multiple edges. 
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2 Description of the algorithm for cubic graphs 

The heuristic algorithm is cahed CVD (for "conflicting vertex displacement") and starts 
with a random coloring of the edges of G using 3 colors. A conflicting vertex is a vertex 
whose incident edges are not properly colored (at least two of them are equally colored). 
The main loop of the algorithm consists of the following steps. First, randomly choose a 
conflicting vertex and take one of the unproperly colored edges as the starting point of a 
Kempe chain. That is, a 2-edge-colored path whose two colors are the original unproper 
color of the starting edge and a new color that makes the vertex unconflicting or, at 
least, reduces the number of edges with repeated color (a suitable new color always exists, 
and can be randomly chosen when it is not unique). Then, swap the two colors along 
the chain. Such an operation does not create new conflicts but moves the conflict along 
the path. Stop the color swapping when, either, another conflicting vertex is reached, 
where the conflicts may cancel out, or the number of swaps attains a given limit (at most 
the number of vertices in the graph). Next, choose again a conflicting (not necessarily 
different) vertex and start a new Kempe chain. Repeat the procedure until, either, all 
conflicts are solved (so getting a 3-edge-coloring) or a fixed limit (made precise below) is 
reached. 

The number of conflicting vertices may stabilize even if G is not a snark (a class two 
graph, see Fiol [2j, Gardner [4J, or Isaacs [6]) because not every color swap succeeds in 
solving one conflict. To avoid an infinite loop, we keep a counter of the number of times 
a new conflicting vertex is chosen and reset this counter to zero only when the number 
of conflicting vertices decreases. If the counter reaches a certain, constant, value R or 
repetition limit then the main loop is restarted from a new initial random coloring of the 
edges. The maximum number of times the main loop can be restarted is the iteration 
limit, L. If we reach this limit, we stop and assume that the cubic graph can not be 
properly 3-edge-colored. Then the non-null probability of making a wrong assumption is 
the price for "breaking" intractability. 

3 Complexity analysis 

Assuming the worst case, every time a conflicting vertex is chosen at most n colors are 
swapped along a Kempe chain, and at most all n vertices in the graph could be conflicting. 
If the number of conflicting vertices ceases to decrease, the choice of a new conflicting 
vertex can be repeated no more than R times, the repetition limit, before the main loop 
restarts from a new random coloring of the edges, and the number of restarts is at most the 
iteration limit, L. The cost of randomly coloring the graph edges is 0{n), and at most L 
random colorings must be done. The parameters R and L are fixed during the execution of 
the algorithm. Therefore, the expected running time is 0{n'^) + 0{n) = 0{n'^). Given the 
pessimistic estimation, it could be better in practice (see Section 6). A similar complexity 
analysis is valid in more general cases (Section 5). 
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4 Algorithm modifications for A-regular graphs and graphs 
with maximum degree A 

The basic idea (conflict displacement) is the same, but this case is slightly more complex. 
We measure the "conflict level" of a vertex (how far it is from having its incident edges 
properly colored) by its degree (A if the graph is A-regular) minus the number of different 
colors used in its incident edges. This gives a measure of the color repetition: assuming 
the colors are always in the range 0,1, 2,..., A — 1, if the number of different colors used 
in the incident edges of a vertex is equal to the degree of that vertex, then no color is 
repeated and the conflict level is 0. If the same color is used in every incident edge, there 
are A — 1 repetitions of that color and the conflict level of the vertex is A — 1. The range of 
possible conflict level values is 0,1, 2,..., A — 1. A data structure, a dictionary to be precise, 
mapping each conflict level value to the set of vertices with that conflict level, is updated 
as the Kempe switches change the coloring of the edges (see Section 8) . This information 
is needed because in the new algorithm conflicting vertices are chosen according to their 
conflict level: vertices with the highest conflict level are chosen first. The conflict level 
value and its corresponding set of unconflicting vertices do not need to be present in the 
dictionary. 

Instead of the number of conflicting vertices, the total number of conflicts is used to 
measure the "conflictivity" of the graph: this is the sum of the conflict levels over all 
vertices in the graph, and it may be larger than the number of vertices in the graph. 

To reduce the conflictivity of the initial arbitrary edge-coloring a simple greedy algorithm 
can be used instead of the random coloring that was sufficient in the cubic case. Besides 
this, and implementation details aside, the conflicting vertex choice based on the highest 
conflict level and the conflictivity measure based on the sum of all conflict levels are the 
only relevant modifications of the algorithm done to deal with the general case. 

5 Complexity analysis in the case of A-regular graphs or 
graphs with maximum degree A 

As in the previous analysis, and asuming the worst case again, every time a conflicting 
vertex is chosen at most n swaps are done along a Kempe chain. In the present case 
we use the sum of all conflict levels to measure the conflictivity of the graph, which is 
at most n(A — 1). If the conflictivity measure ceases to decrease, the choice of a new 
conflicting vertex can be repeated no more than R times, the repetition limit, before the 
main loop restarts from a new pre-coloring of the edges, and the number of restarts is at 
most the iteration limit, L. The cost of pre-coloring all 0(An) graph edges using a random 
coloring or the greedy algorithm is O(A^n) or possibly better, depending on the efficiency 
of sampling from the set of 0(A) available colors, and at most L pre-colorings must be 
done. The parameters R and L are fixed during the execution of the algorithm. Therefore, 
the expected running time is 0((A — l)n^) + O(A^n) = O(A^n^). As in the previous case. 
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the estimation is pessimistic and the performance could be better in practice (see Section 
7). 

Roughly speaking, our approach, with its vertex-centered measure of conflict, is similar 
to the edge-centered approach by Lee, Wan and Guan ([7]). These authors improved on 
previous A-edge-coloring algorithms by giving a polynomial (albeit probabilistic) algo- 
rithm. The CVD algorithm is simpler and easy to implement (see Section 8). 



6 Performance testing in the case of cubic graphs 

We tested the performance of the heuristic algorithm by measuring the time needed to 
3-edge-color random cubic graphs. More specifically, for each randomly generated cubic 
graph, we measured the time needed to find a 3-edge-coloring, the number of iterations 
of the main loop, and the average time per iteration (the coloring time divided by the 
number of iterations). For each number of vertices, we obtained the minimum, average, 
and maximum values of each of these magnitudes over a set of 30 random instances, and 
we plotted the results. It seems that the heuristic algorithm finds a 3-edge-coloring, if any, 
with high probability. 
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The average running time appears to grow linearly with the number of vertices (see the 
first figure). The same happens with the average time per iteration (not shown), whereas 
the average number of iterations does not seem to depend on the number of vertices (see 
the second figure). A linearly growing average running time is in agreement with the 
O(n^) expected running time. 



7 Performance testing in the case of A-regular graphs 

We followed the same testing procedure as in the cubic case (measuring the average running 
time over 30 instances), but varying the number n of vertices and the value of A. The 
growth of the average running time is in agreement with the expected O(A^n^) behaviour. 
The following figure shows four plots, corresponding to four different values of the variable 
A (3, 7, 11 and 15). 
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8 Implementation details 

We used the Python language (version 2.7.3) to implement the heuristic algorithm. This 
language allows us to use the NetworkX package to generate random regular graphs, and 
the matplotlib package to plot the results of the performance tests. 

In order to facilitate the reproducibility of our results, or the use of the CVD algorithm to 
anyone interested, we give a complete Python implementation of the algorithm described 
in Section 4, depending only on the Python Standard Library and the NetworkX package. 
The plot of Section 6 was created using the algorithm described in Section 3, specific for 
cubic graphs, with a very similar (slightly simpler) implementation. 

### Implementation data ### 

#Operating system: Linux-3 . 2 . 0-31-generic-x86_64-with-Ubuntu-12 . 04-precise 
#CPU: Intel (R) Core(TM) 13 CPU M 370 2.40GHz. Cache size: 3072 KB 

#CPU speed: 933.000 Hz 

### Python code ### 

from random import choice, randint 

import networkx as nx 
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def properlyColored(G,u,D) : 

return lenCset (G [u] . values ()))==G. degree (u) and alKcolor in range(D) for 

color in G[u] .values ()) 

def checkEdgeColoring(G,D) : 

return all(properlyColored(G,u,D) for u in G.nodesO) 

def conf lictLeveKGjU) : return G.degree(u)-len(set(G[u] .valuesO)) 

def createConf lictDictionary(G,D) : 

conflict_dictionary=dict( [(i,set( □ )) for i in range (1,D)]) 
for u in G.nodesO: 

conf lict_level_u=conf lictLevel (G ,u) 

if conf lict_level_u>0 : conf lict_dictionary [conf lict_level_u] .add(u) 
return conf lict_dictionary 

def updateConf lictDictionary (G,u,conf lict_dictionary ,old_conf lict_level_u) : 
conf lict_level_u=conf lictLevel (G,u) 
if old_conf lict_level_u>0 : 

conf lict_dictionary [old_conf lict_level_u] . remove (u) 
if conf lict_level_u>0: 

conf lict_dictionary [conf lict_level_u] . add(u) 
return conf lict_level_u-old_conf lict_level_u 

def maxConf lictLevel (conf lict_dictionary) : 

return max( [conf lict_level for conf lict_level in conf lict_dictionary 
if len(conf lict_dictionary [conf lict_level] )>0] ) 

def totalNumberOf Conf licts (conf lict_dictionary) : 

return sum(conf lict_level*len(conf lict_dictionary [conf lict_level] ) 
for conf lict_level in conf lict_dictionary) 

def colorEdgeAndUpdate(G,u,v, color, conf lict_dictionary) : 
old_conflict_level_u=conf lictLevel (G,u) 
old_conflict_level_v=conf lictLevel (G,v) 
G[u] [v]=G[v] [u]=color 
updateConf lictDictionary (G,u, 

conf lict_dictionary , 
old_conf lict_level_u) 
return updateConf lictDictionary (G , v , 

conf lict_dictionary, 
old_conf lict_level_v) 
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def KempeNext(G, last, node, new_color,conflict_dictionary) : 
available_f or_next= [w for w in G[node] if w!=last 

and G [node] [w] ==new_color] 
if available_f or_next== [] : next_node=None 
else : next_node=choice (available_f or_next) 
old_color=G [last] [node] 

conf lict_level_variat ion=colorEdgeAndUpdate (G , last , node , new_color , 

conflict .dictionary) 
return conf lict_level_variation,old_color , next _node 

def KempeStepCG, last, node, new_color, conf lict_dictionary) : 

conf lict_level_variation, old_color , next _node=KempeNext (G , last , node , 

new_color , 

conf lict_dictionary) 
if conf Iict_level_variation<0 or next_node==None : return node , None , None 
return node, next _node,old_color 

def KempeProcess(G, last, node, new_color, conf lict_dictionary) : 
Kempe_chain=set ( [] ) 

while new_color!=None and last not in Kempe_chain: 
Kempe_chain. add(last) 

last , node , new_color=KempeStep (G , last , node , new_color , 

conf lict_dictionary) 

def KempeStart(G,D,node,conflict_dictionary) : 
colors=set (range (D) ) 
next_node=None 
for adjacent in G[node]: 

edge_color=G [node] [adjacent] 

if edge_color in colors: colors. remove (edge_color) 
else: next_node=adjacent 
if next_node ! =None : 

KempeProcess (G , node , next _node , choice (list (colors) ) , conf lict_dictionary) 

def preColoring(G,D) : #Pre-coloring with a greedy algorithm 
for e in G.edgesO: G[e[0]] [e[l]]=G[e[l]] [e[0]]=None 
for e in G.edgesO: 

available_colors=set (range (D) ) 

available_colors-=set (G [e [0] ] . values () ) 

available_colors-=set (G [e [1] ] . values () ) 

if available_colors==set() : G[e[0]] [e [1] ] =G [e [1] ] [e [0]]=randint(0,D-l) 
else : G [e [0] ] [e [1] ] =G [e [1] ] [e [0] ] =choice (list (available_colors) ) 
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#def preColoring(G,D) : #Random pre-coloring 

# for e in G.edgesO: G [e [0] ] [e [1] ] =G [e [1] ] [e [0] ] =randiiit (0 ,D-1) 

def heuristic(G,D,repetition_liinit) : 
repetitioiiCounter=0 

conf lict_dictionary=createConf lictDictionary (G,D) 
previous=curr ent=totalNumber Of Conflict s( conflict _dictionary) 
while previous>0: 

highest _ conf lict_level=niaxConflictLevel (conf lict_dictionary) 
node=choice (list (conf lict_dictionary [highest_conf lict_level] ) ) 
KempeStart (G, D, node, conf lict_dictionary) 
current=totalNumberOf Conf licts (conf lict_dictionary) 
if current==0 : return True 
if current>=previous : 

repetitionCounter+=l 

if repetitionCounter>repetition_liiiiit : return False 
else: repetitionCounter=0 
previous=min (previous , current) 
return True 

def applyHeuristic(G,D,repetition_limit , iteration_limit) : 
preColoring (G , D) 
number_of _iterations=l 

while not heuristic(G,D,repetition_limit) : 

if number_of _iterations>iteration_limit : break 

preColoring (G , D) 

number _of_iter at ions+=l 
print "Number of iterations :" ,number_of_iterations 
print "Edge-coloring successful: " ,checkEdgeColoring(G,D) 

### Example ### 

repetition_limit=iteration_limit=50 

NUMBER_OF_VERTICES=50000 

DEGREE=3 

G=nx.random_regular_graph (DEGREE, NUMBER_OF_VERTICES) 
applyHeuristic (G, DEGREE, repet it ion_limit , iteration_limit) 



9 Testing a conjecture on odd graphs 

Odd graphs, 0^ with A: > 2 an integer, are defined in the following way: the vertices 
correspond to the {k — l)-subsets of a {2k — l)-set, and two vertices are adjacent if their 
corresponding subsets are disjoint. It is conjectured (Biggs Fiorini and Wilson [3]) 
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that odd graphs are of class one except for k = 3 (the Petersen graph) and k a power of 
two (when the graph has an odd number of vertices, implying that it is trivially of class 
two). In their book, Fiorini and Wilson comment results for k < 8. We used the CVD 
algorithm to test that the conjecture is true for A; < 11. 

10 Conclusions and future directions 

We have seen that random cubic and A-regular graphs can be edge-colored with empirical 
efficiency by a relatively simple heuristic algorithm (with simple data structures). 

Many combinatorial problems, for instance logical circuit problems (Fiol [2]), can be 
reduced to 3-edge-coloring of cubic graphs or edge-coloring of regular graphs. Therefore, 
any efficient edge-coloring heuristic might be useful to solve them. Moreover, efficient 
edge-coloring is of practical interest in many applications. 

We are currently working in the generalization of the heuristic algorithm to hypergraphs 
and other combinatorial structures, and testing its performance in combinatorial problems. 
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